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Ankudinov  Ship  Squat  Predictions  - 
Part  I:  Theory,  Parameters,  and 
us  Army  corps  FORTRAN  Programs 

of  Engineers® 

by  Michael  J.  Briggs 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  summarizes 
the  Ankudinov  empirical  formula  for  ship  squat  predictions  used  by  the  U.S.  Army  Engineer 
Research  and  Development  Center  (ERDC)  Ship  Tow  Simulator  (STS).  This  CHETN  documents 
the  Ankudinov  squat  formulas  and  describes  two  FORTRAN  programs  that  were  written  for  sin¬ 
gle  and  multiple  ship  speed  applications.  Two  examples  are  given  for  a  Panamax  bulk  carrier  and 
Panamax  tanker  for  both  programs.  A  companion  CHETN  compares  and  validates  the  Ankudi¬ 
nov  squat  predictions  with  laboratory  measurements  of  a  Post-Panamax  containership  and  field 
measurements  of  a  Panamax  containership,  tanker,  and  bulk  carrier  in  the  Panama  Canal.  It  also 
compares  the  Ankudinov  formula  with  several  of  the  Permanent  International  Association  of 
Navigation  Congresses  (PIANC)  empirical  squat  formulas. 

BACKGROUND:  Squat  is  the  reduction  in  underkeel  clearance  (UKC)  between  a  vessel  at-rest 
and  underway  due  to  the  increased  flow  of  water  past  the  moving  hull.  The  forward  motion  of 
the  ship  pushes  water  ahead  of  it  that  must  return  around  the  sides  and  under  the  keel.  This  water 
motion  induces  a  relative  velocity  between  the  ship  and  the  surrounding  water  that  causes  a 
water  level  depression  in  which  the  ship  sinks.  The  velocity  field  produces  a  hydrodynamic  pres¬ 
sure  change  along  the  ship  similar  to  the  Bernoulli  effect.  This  phenomenon  produces  a  down¬ 
ward  vertical  force  (sinkage,  positive  downward)  and  a  moment  about  the  transverse  axis  (trim, 
positive  bow  up)  that  can  result  in  different  values  of  squat  at  the  bow  and  stern  (Figure  1).  This 
combination  of  sinkage  and  change  in  trim  is  called  ship  squat. 

The  main  ship  parameters  include: 

•  Ship  draft  or  draught  T 

•  Hull  shape  as  represented  by  the  block  coefficient  Cb 

•  Ship  speed  Vs  (m/sec)  or  Vk  (knots) 


Figure  1 .  Schematic  of  ship  squat  at  bow  and  stern  (figure  courtesy  of  Handbook  of  Coastal  and  Ocean 

Engineering,  World  Scientific  Publishers). 
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Other  ship  parameters  include  the  length  between  forward  and  aft  perpendiculars  Lpp  and  the 
beam  B.  The  Cb  is  a  measure  of  the  “fineness”  of  the  vessel’s  shape  relative  to  an  equivalent 
rectangular  volume  with  the  same  dimensions.  The  range  of  values  of  Cb  is  typically  between 
0.45  for  high-speed  vessels  and  0.85  for  slow,  full-size  tankers  and  bulk  carriers.  Typically,  full- 
form  ships  like  tankers  and  bulk  carriers  squat  by  (i.e.,  at)  the  bow  Sb  and  more  slender  ships  like 
containerships  by  the  stem  Ss.  The  initial  trim  of  the  ship  also  influences  the  location  of  the 
maximum  squat.  The  most  important  ship  parameter  is  its  speed  Vs-  The  main  channel  consid¬ 
erations  are  proximity  of  the  channel  sides  and  bottom  and  cross-sectional  configuration.  The 
three  channel  configurations  are  unrestricted  or  open  or  fairway  (U),  restricted  or  trench  (R),  and 
canal  (C).  Channel  parameters  include  channel  depth  h,  inverse  side  slope  n  (i.e., 
n  =  run/rise  =  cot  9),  and  trench  height  from  the  bottom  of  the  channel  to  the  top  of  the  trench  hr- 
For  a  canal  cross  section,  hr  equals  h. 

This  CHETN  builds  on  the  earlier  CHETN  entitled  “Ship  Squat  Predictions  for  Ship/Tow  Simu¬ 
lator,”  CHETN-I-72  (Briggs  2006)  that  summarized  several  PIANC  empirical  ship  squat  predic¬ 
tions.  The  Ankudinov  squat  predictions  are  currently  used  in  the  Corps  STS.  In  the  first  section, 
the  equations  are  described  for  the  Ankudinov  ship  squat  predictions.  In  the  next  section,  the 
FORTRAN  programs  ANKUDINOV4  for  single  ship  speed  and  ANKUDINOVM4  for  multiple 
ship  speeds  are  described.  An  application  of  CHETN-I-63  (Demirbilek  and  Sargent  1999)  and 
CHETN-I  -72  (Briggs  2006),  the  Panamax  bulk  carrier  is  used  to  describe  the  input  and  output 
parameters  and  results  from  the  Ankudinov  predictions  relative  to  PIANC  (1997)  empirical  for¬ 
mulas  for  unrestricted  (open)  applications.  Finally,  in  Engineer  Manual  (EM)  1110-2-1613  the 
Panamax  tanker  is  used  as  an  example  to  illustrate  the  effect  of  channel  type  on  the  Ankudinov 
predictions  for  bow  squat.  Part  II  of  this  CHETN  will  compare  and  validate  the  Ankudinov  pre¬ 
dictions  with  laboratory  and  field  measurements  for  several  ship  and  channel  combinations. 

ANKUDINOV  SQUAT  FORMULA:  Ankudinov  and  Jakobsen  (1996)  and  Ankudinov  et  al. 
(1996,  2000)  proposed  the  MARSIM  2000  formula  for  maximum  squat  based  on  a  midpoint 
sinkage  Sm  (SMid)  and  vessel  trim  Trim  (Trim)  in  shallow  water.  The  names  in  parenthesis  cor¬ 
respond  with  the  nomenclature  used  in  the  FORTRAN  programs.  Appendix  A  contains  a  list  of 
all  the  parameters  in  the  Ankudinov  programs  for  ease  of  reference.  The  Ankudinov  method  has 
undergone  considerable  revision  as  new  data  were  collected  and  compared.  The  most  recent 
modifications  from  a  study  of  ship  squat  in  the  St.  Lawrence  Seaway  (Stocks  et  al.  2002),  are 
e-mails  and  telecons  in  April  2009  (Ankudinov  2009)  that  are  contained  in  this  CHETN  and  the 
FORTRAN  programs. 

The  Ankudinov  prediction  is  one  of  the  most  thorough,  but  also  the  most  complicated  formulas 
for  predicting  ship  squat.  These  components  include  factors  to  account  for  the  effects  of  the  ship 
and  channel.  The  restriction  on  Depth  Froude  Number  Fnh  are  values  less  than  or  equal  to  0.6. 
The  maximum  ship  squat  Smox  (SMax)  is  a  function  of  two  main  components:  the  midpoint  sink- 
age  Sm  (SMid)  and  the  vessel  trim,  Trim  (Trim)  given  by: 

S^=LjS.  +  0.5Trim)  (1) 

The  Smcix  can  be  at  the  bow  or  stern  depending  on  the  value  of  Trim.  The  negative  sign  in  Equa¬ 
tion  1  is  used  for  bow  squat  Sb  (SbA)  and  the  positive  sign  for  stern  squat  Ss  (SsA). 
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Midpoint  Sinkage  Sm  (SMid).  The  Sm  (SMid)  is  defined  as: 

S„=SMid  =  (l  +  Ksp)PH„PFPthlTPcu  (2) 

The  ship,  water  depth,  and  channel  parameters  in  this  midpoint  sinkage  equation  are  described  in 
the  following  paragraphs.  The  propeller  parameter  Ksp  (KpS)  is  defined  as: 

s  f0. 1 5  single  propeller 

Kp  =  KpS  =  <  (3) 

[0.13  twin  propellers 

The  ship  hull  parameter  for  shallow  water  PHu  (Par_Hull_s)  was  recently  modified  by  Ankudi¬ 
nov  (2009)  as: 


PHu  =  Par  _  Hull  _s  =  \JCB 


f  bt^ 

y^ppj 


+  0.004C„ 


The  ship  forward  speed  parameter  Ppnh  (Par_Fnh)  is  given  by: 


Pr,  =  Par  Fnh  =  F, 


_  „(1.8+0.4F„t) 


nh 


(4) 


(5) 


which  is  a  numerical  approximation  to  the  tenn  “  FiBi /yjl-  F„h  ”  that  is  in  many  of  the  PIANC 
empirical  squat  fonnulas.  The  Fnh  (Fnh)  is  defined  as: 


Fnh  =  Fnk  = 


V. 


sfgh 


(6) 


The  water  depth  effects  parameter  P+h/r  (Par_hT_s)  is  defined  as: 


P+h/T  =  Far  _  hT  _  s 


1.0  + 


0.35 

(hlrf 


(7) 


Note  that  the  h/T  ratio  corresponds  to  the  nondimensional  ratio  R/,r- 
The  channel  effects  parameter  Pchi  (P_Chl)  is  given  by: 


PChl=P_Chl  = 


[1.0  u 

[L0  +  l0Sh-l.5(l.0  +  Sh)^  R,C 


(8) 


where  the  canal  or  restricted  configuration  is  incorporated  in  the  channel  depth  factor  S/,  (Sh) 
defined  by: 
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f  c  V(,  \ 

S,=Sh  =  CB  —  R,C  (9) 

B{h/Tj{h) 

and  the  blockage  factor  S  (=  AfAch )  is  the  fraction  of  the  cross-sectional  area  of  the  waterway 
A  a,  that  is  occupied  by  the  ship’s  underwater  midships  cross  section  As. 

Vessel  Trim.  The  second  main  component  in  the  MARSIM  squat  equation  is  the  vessel  trim. 
Trim  that  was  also  recently  modified  by  Ankudinov  (2009)  as: 

Trim  =  -\.lPHuPFi/hjTKTrPChl  (10) 

Note  that  the  LPPIB  ratio  corresponds  to  the  nondimensional  ratio  Rlb-  In  addition  to  the  three 
parameters  already  described  for  the  midpoint  sinkage  equation,  the  Trim  also  includes  parame¬ 
ter  Ph/T  (ParhT)  and  coefficient  KTr  (KTr)  to  quantify  the  effects  of  the  ship  propellers,  bulbous 
bow,  stern  transom,  and  initial  trim. 

The  vessel  trim  parameter  Pyr  (Par  hT)  accounts  for  the  reduction  in  trim  due  to  the  propeller  in 
shallow  water  and  is  defined  as: 

"2.5(1  -h/T)' 

Ph/T  =  Par  _hT  =  l-e'~  F"h  *  (11) 

The  Trim  coefficient  Kjr  (K_Tr)  is  a  function  of  many  factors  and  is  given  by: 

KTr  =  K _Tr  =  Cf  - (0. 1 5KSP  +  KTP  )  - (ktb  +  KTTr  +  Krn  )  (12) 

The  first  factor  in  this  equation  C’f  '  is  the  block  coefficient  Cg,  raised  to  the  njr  (nTr)  power. 
This  Trim  exponent  nr,-  is  defined  as: 

nTr=nTr  =  2.0  + 0.8^-  (13) 

CB 

The  next  two  factors  define  the  propeller  effect  on  the  vessel  trim.  The  first  factor  Ksp  (KpS)  is 
the  same  as  the  propeller  parameter  for  the  midpoint  sinkage  and  the  second  factor  is  the  pro¬ 
peller  trim  parameter  KTB  (KpT): 

T  f 0. 1 5  single  propeller 

KTp  =  KpT  =  \  Bee  (14) 

[0.20  twin  propellers 

The  last  group  of  three  factors  define  the  effects  of  the  bulbous  bow  K[  (KbT),  stern  transom 
KTTr  (KTrT),  and  initial  trim  K.'n  (KT1T)  on  the  vessel  trim.  The  bulbous  bow  factor  KTb  (KbT) 
is  given  by: 
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K[  =  KbT  = 


0.1 


0.0 


bulbous  bow 
no  bulbous  bow 


The  Stern  transom  factor  KTTr  (KTrT)  is  defined  by: 


(15) 


0.1 

Bn 

=  0.1 

0.45 

KTTr  =  KTrT  =  < 

B  J 

B 

0.0 

stem  transom 
no  stern  transom 


(16) 


where  BTr  (BTr)  is  the  Stern  transom  width  and  is  typically  0.4  B,  although  values  as  high  as 
0.7  B  have  sometimes  been  used. 

The  initial  trim  effect  factor  K7n  (KT1T)  is  given  by: 


K'n  =  KTIT  = 


(A+A) 


(17) 


where  Tap  (Tap)  is  the  static  draft  at  the  stern  or  aft  perpendicular  and  TfP  (Tfp)  is  the  static  draft 
at  the  bow  or  forward  perpendicular. 

Finally,  the  channel  effect  trim  correction  parameter  (Par_Ch2)  Paa  is  defined  as: 


p  =  Par 

1  Chi  1  ur 


Ch2  = 


P 

[1-0-5S, 


u 

R,C 


(18) 


FORTRAN  PROGRAMS  ANKUDINOV4  and  ANKUDINOVM4 :  The  two  FORTRAN  pro¬ 
grams  ANKUDINOV4  and  ANKUDINOVM4  are  similar  to  each  other  as  the  one  with  the  “M” 
in  the  name  (for  multiple  ship  speeds)  is  a  copy  of  the  other  for  single  ship  speeds.  They  both 
reflect  the  latest  Version  3  modifications.  The  multiple  speed  version  can  have  any  beginning, 
increment,  and  number  of  speeds  as  long  as  not  greater  than  100.  In  2000,  both  programs  were 
written  in  Compaq  Visual  Fortran,  Version  6.5.  They  are  similar  to  the  program  SQUAT  for  the 
PIANC  formulas  described  in  CHETN-I-72  (Briggs  2006).  They  are  modular  as  most  compo¬ 
nents  are  contained  in  separate  subroutines.  The  user  can  enter  input  data  in  two  different  ways: 
in  manual  entry  mode  (Subroutine  Input)  or  by  reading  (Subroutine  Readln)  a  previously  stored 
file  of  the  input  variables.  Both  are  perfonned  in  real  time,  but  the  Readln  option  is  more  effi¬ 
cient  as  the  user  does  not  have  to  re-enter  input  for  each  ran.  Input  data  are  automatically  stored 
(Subroutine  Store)  in  a  generic  input  file  Ankudln.out  (AnkudMIn.out  for  ANKUDINOVM4) 
after  input  is  completed.  The  output  is  stored  in  a  generic  file  named  AnkudOut.out  (Ankud- 
MOut.out  for  ANKUDINOVM4).  The  user  can  (and  should)  rename  these  generic  input  and 
output  files  to  save  the  data  for  later  use  and  analysis. 

The  program  is  well  documented  with  comments  throughout.  The  naming  convention  follows  the 
PIANC  (1997)  and  Ankudinov  et  al.  (1996,  2000)  nomenclature.  Units  are  in  the  metric  system 
with  predicted  squat  in  meters.  The  program  is  completely  stand-alone  at  this  point  so  that  it  can 
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be  used  in  comparisons  with  other  squat  predictions  (see  CHETN-I-72).  Interested  readers  can 
obtain  a  copy  of  the  executable  version  of  the  program  and  I/O  files  by  contacting  the  author. 

Input.  The  first  input  question  is  whether  the  user  wants  to  manually  input  (enter  I)  or  read  in 
(enter  R)  the  data  from  a  previously  stored  file.  The  next  question  is  the  type  of  channel.  The 
user  can  enter  U  for  unrestricted,  R  for  restricted,  or  C  for  canal.  The  input  is  the  same  for  all 
three  channel  configurations  for  most  of  the  questions,  but  each  has  some  slight  variations.  These 
first  two  inputs  are  all  in  capital  letters.  If  the  user  accidentally  enters  in  small  letters,  the  pro¬ 
gram  will  abort.  The  user  must  then  restart  (CNTL  C  keys  together  to  cancel  the  open  DOS  win¬ 
dow)  and  re-enter  the  input  values  in  capital  letters. 

The  next  questions  involve  the  ship  and  channel  parameters.  Table  1  lists  input  parameters,  sym¬ 
bols,  and  units  for  both  programs  for  all  three  channel  configurations.  Figure  2  is  an  example  for 
an  unrestricted  channel  (U)  for  (a)  ANKUDINOV4  and  (b)  ANKUDINOVM4  programs  using 
the  R  (i.e.,  Readln  option).  This  example  is  for  the  modified  Bunga  Saga  Emphat  Panamax  bulk 
carrier  that  was  originally  described  by  Harkins  and  Dorrell  (2000)  and  used  in  CETN-I-63  and 
CHETN-I-72.  Table  2  summarizes  the  input  parameters  for  this  example  as  well  as  the  R  and  C 
channel  types.  After  completing  input,  the  program  runs  without  any  further  action  from  the  user. 


Table  1.  ANKUDINOV4  and  ANKUDINOVM4  input  parameters. 

Parameter  Description 

Symbol 

Units 

Program 

Channel  Type 

Unrestricted  U 

Restricted  R 

Canal  C 

Length  between  perpendiculars 

Lpp 

m 

Both 

Y 

Y 

Y 

Beam 

B 

m 

Both 

Y 

Y 

Y 

Draft 

T 

m 

Both 

Y 

Y 

Y 

Draft  at  forward  perpendicular 

TfP 

m 

Both 

Y 

Y 

Y 

Draft  at  aft  perpendicular 

Tap 

m 

Both 

Y 

Y 

Y 

Block  coefficient 

CB 

— 

Both 

Y 

Y 

Y 

Vessel  speed  or  1st  speed 

vk 

knots 

ANKM 

Y 

Y 

Y 

Vessel  increment  speed 

Vmc 

knots 

ANKM 

Y 

Y 

Y 

Number  of  speeds 

NVel 

— 

ANKM 

Y 

Y 

Y 

Number  propeller  flag 

IProp 

— 

Both 

Y 

Y 

Y 

Bulbous  bow  flag 

IBow 

— 

Both 

Y 

Y 

Y 

Transom  stern  flag 

IStern 

— 

Both 

Y 

Y 

Y 

Channel  width  at  bottom 

W 

m 

Both 

N 

Y 

Y 

Bank  slope  (inverse) 

n 

— 

Both 

N 

Y 

Y 

Height  dredged  underwater  trench 

hT 

m 

Both 

N 

Y 

N 

Notes: 

1 .  Both  =  Both  Programs  ANKUDINOV4  and  ANKUDINOVM4. 

2.  AN  KM  =  Program  ANKUDINOVM4  only. 

3.  N  =  No,  Y  =  Yes. 
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251.1600 
32.2500 

12.8000 

12.8000 

12.8000 
0.9050 
10.0000 
15.  3600 
1 
1 
0 

panamax  Bulk  carrier,  cetn-i-63,  Mar  99,  Demirbilek  &  Sargent 
unrestricted  channel,  Ankudinov4  -  29Apr09  changes,  29Apr09 
|  25171600 

32.2500 

12.8000 

12.8000 

12.8000 
0.9050 
5.0000 
1.0000 
11.0000 
15.  3600 
1 
1 
0 

Panamax  Bulk  carrier,  CETN-I-63,  Mar  99,  Demirbilek  &  Sargent 
Unrestricted  channel,  AnkudinovM4  -  29Apr09  changes,  29Apr09 


vessel  length  between  perpendiculars  Lpp,  m 

vessel  beam  B,  m 

vessel  draft  t,  m 

vessel  draft  TFp,  m 

vessel  draft  TAp,  m 

vessel  block  coefficient  cb 

vessel  speed  vk,  kts 

water  depth  h,  m 

single  (1)  or  twin  (2)  propeller,  iProp 
Bulbous  bow  (1)  or  conventional  (0),  ibow 
Stern  transom  (1)  or  conventional  (0),  IStern 


vessel  length  between  perpendiculars  Lpp.  m 

vessel  beam  b,  m 

vessel  draft  T,  m 

vessel  draft  TFp,  m 

vessel  draft  TAp,  m 

vessel  block  coefficient  CB 

vessel  speed  vk,  kts 

vessel  increment  speed  vine,  kts 

Number  vessel  speeds  NVel 

water  depth  h,  m 

single  (1)  or  twin  (2)  propeller,  IProp 
Bulbous  bow  (1)  or  conventional  (0),  ibow 
stern  transom  (1)  or  conventional  (0),  IStern 


Figure  2.  Example  generic  input  files  for  modified  Bunga  Saga  Emphat  Panamax  bulk  carrier  and  unre¬ 
stricted  channel  for  (a)  single  ship  speed  Ankudln.out  and  (b)  multiple  ship  speeds  AnkudMIn.out. 


Table  2.  Ship  and  channel  input  parameters  for  modified  Bunga  Saga  Emphat 

Panamax  bulk  carrier. 

Lpp  (m) 

e(m) 

T  (m) 

CB 

h(m) 

M  m) 

W  (m) 

n 

251.2 

32.3 

12.8 

0.91 

15.4 

3.84 

281 

3 

Notes: 

1.  Tfp  =  Tap  =  T;  IProp  =  1;  IBow  =  1;  IStern  =  0;  KpS 

=  0.15;  KpT  =  0.15. 

Output.  The  output  is  first  printed  to  the  screen  in  a  DOS  Window  for  the  user  to  check.  The 
user  can  examine  the  contents  as  long  as  necessary.  A  carriage  return  will  close  this  window.  The 
output  is  then  automatically  written  to  the  generic  output  file  AnkudOut.out  (AnkudMOut.out  for 
ANKUDIN  O  VM4) . 

Figures  3a  and  b  show  the  output  corresponding  to  the  inputs  in  Figure  2  for  ANKUDINOV4 
and  ANKUDINOVM4,  respectively.  The  project  title,  ship  input  parameters,  channel  input 
parameters,  and  calculated  parameters  are  identical  for  both  programs.  The  output  parameters 
consisting  of  ratios  and  constants  and  the  Ankudinov  output  parameters  that  include  K  factors 
and  constants,  Par  factors  and  constants,  and  sinkage  and  trim  results  are  next  printed.  The  main 
difference  between  the  programs  is  that  the  parameters  that  are  affected  by  ship  speed  are  listed 
separately  as  a  function  of  speed.  Also,  the  njr  trim  exponent  is  not  included  for  the  multiple 
speed  version. 

Results.  Figure  4  is  a  plot  of  the  Ankudinov  bow  Si,  (SbA)  and  stem  Ss  (SsA)  squat  predictions 
for  the  modified  Bunga  Saga  Emphat  bulk  carrier  in  an  unrestricted  channel.  Averages  from  the 
PIANC  empirical  fonnulas  (CHETN-I-72)  are  also  shown  for  reference.  Table  3  lists  the  pre¬ 
dicted  values  in  meters  for  the  Ankudinov  programs  at  5  and  10  knots  for  all  three  channel  types 
U,  R,  and  C.  The  non-dimensional  ratio  of  the  Ankudinov  to  PIANC  predictions  at  bow  and  stem 
are  also  included  for  comparisons. 
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Results  From  Program  Ankudinov4,  version  4.0 

panamax  Bulk  carrier,  cetn-i-63,  Mar  99,  Demirbilek 
unrestricted  channel,  Ankudinov4  -  29Apr09  changes, 

&  sargent 

29Apr09 

user  accepts  waiver  of  warranty,  Limitation  of  Liability,  indemnity,  and  Assent 

contained  in  the  ReadMe  file  for  program  Ankudinov4 

and  its  versions 

***  ship  input  Parameters  *** 

Length  between  perpendiculars,  Lpp  = 

251.16  m 

Beam,  B  = 

Draft,  T  = 

32.25  m 

12.80  m 

Draft  Forward  perpendicular.  Tfp  = 

12.80  m 

Draft  Aft  Perpendicular,  Tap  = 

12.80  m 

Block  coefficient,  CB  = 

0.905 

***  channel  input  Parameters  *** 

channel  type  code  CHType  = 
c  =  canal 

R  =  Restricted  channel  w/trench 
u  =  unrestricted  channel  or  flat 

U 

water  depth,  h  = 

15.36  m 

channel  effective  width,  WEff3  = 
used  for  channel  width 

247.14  m 

***  calculated  Parameters  *** 

ship  displacement  volume,  vol  =  93829.36  mM 

coef  for  as  bilge  radius,  cs  = 

0.98 

Area  ship  midship,  As  = 

Area  channel  or  canal,  ACh  = 

404. 54  m^2 

3796.14  mA2 

Blockage  factor,  s  = 

***  output  Parameters  *** 

Ratios  and  constraints 

0.107 

Block  coefficient,  CB  = 

0.905 

ship  speed,  vk  = 

10.00  kts 

ship  speed,  vs  = 

5.14  m/s 

Depth  Froude  No. ,  Fnh  = 

Ratio  depth  to  draft,  RhT=h/T  = 

0.419 

1.20 

Ratio  ship  length  to  depth,  RLh=Lpp/h  = 

16.  35 

Ratio  ship  length  to  beam,  RLB=Lpp/B  = 

7.79 

Ratio  ship  length  to  draft,  RLT=Lpp/T  = 

19.  62 

Ratio  ship  beam  to  draft,  rbt=b/t  = 

***  Ankudinov  output  Parameters  *** 

K  Factors  and  constants 

2.  52 

Single/Twin  Propeller  sinkage,  Kps  = 

0.15 

channel  factor,  sh  = 

0.080 

channel  effect.  P  chi  = 

1.000 

channel  trim  effect,  P_ch2  = 

1.000 

Trim  exponent,  nTr  = 

2.88 

Single/Twin  propeller  trim,  KpT  = 

0.15 

Bulbous  Bow  trim,  KbT  = 

0.10 

stern  transom  trim,  KTrT  = 

0.00 

initial  trim,  KTlt  = 

0.00 

Trim  coefficient,  K_Tr  = 

0.477 

Par  Factors  and  constants 

Hull  effect  shallow,  Par  Hull  s  = 

0.0133 

water  depth  effect,  Par  hT  s  = 

1.2431 

Forward  speed  effect,  par  Fnh  = 

0.1806 

Propeller  effect  shallow,  Par_hT  = 

0. 6967 

sinkage  &  Trim  Results 

Mid-snip  sinkage,  SMid  = 

0.00345  m/m 

ship  Trim,  Trim  = 

0.00136  m/m 

Bow  squat,  sbA  = 

1.04  m 

stern  squat,  ssA  = 

0.  69  m 

Figure  3a.  Example  generic  output  files  for  modified  Bunga  Saga  Emphat  Panamax  bulk  carrier 
and  unrestricted  channel  for  single  ship  speed  AnkudOut.out. 
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Results  From  program  AnkudinovM4,  version  4.0 

Panamax  Bulk  carrier,  cetn-i-63,  Mar  99,  Demirbilek  &  Sargent 

unrestricted  channel,  AnkudinovM4  -  29Apr09  changes,  29Apr09 

user  accepts  waiver  of  warranty,  Limitation  of  Liability,  indemnity,  and  Assent 
contained  in  the  ReadMe  file  for  Program  AnkudinovM4  and  its  versions 


***  ship  input  parameters  *** 
Length  between  perpendiculars,  Lpp 
Beam,  B  = 

Draft,  T  = 

Draft  Forward  Perpendicular,  Tfp  = 
Draft  Aft  Perpendicular,  Tap  = 
Block  coefficient,  CB  = 

***  channel  input  Parameters  *** 
channel  type  code  CHType  = 
c  =  canal 

R  =  Restricted  channel  w/trench 
u  =  unrestricted  channel  or  flat 
water  depth,  h  = 
channel  effective  width,  WEff3  = 
used  for  channel  width 


251.16  m 
32.25  m 
12.80  m 
12.80  m 
12.80  m 
0.905 


u 


15.36  m 
247.14  m 


***  Calculated  Parameters  *** 
ship  displacement  volume,  vol  = 
coef  for  As  bilge  radius,  cs  = 
Area  ship  midship,  as  = 

Area  channel  or  canal,  ACh  = 
Blockage  factor,  s  = 


93829. 36  mA3 
0.98 

404. 54  m*2 
3796.14  mA2 
0.107 


***  output  Parameters  *** 

Ratios  and  constraints 

Ratio  depth  to  draft,  RhT=h/T  = 

Ratio  ship  length  to  depth,  RLh=Lpp/h  = 
Ratio  ship  length  to  beam,  RLB=Lpp/B  = 
Ratio  ship  length  to  draft,  RLT=Lpp/T  = 
Ratio  ship  beam  to  draft,  rbt=b/t  = 


1.  20 
16.  35 
7.79 
19.62 
2.  52 


***  Ankudinov  output  Parameters  *** 

K  Factors  and  constants 

Single/Twin  Propeller  sinkage,  Kps  =  0.15 

channel  factor  effect,  sh  =  0.080 

channel  effect,  P_chl  =  1.000 

channel  trim  effect,  P_ch2  =  1.000 

Single/Twin  Propeller  trim,  KpT  =  0.15 

Bulbous  bow  trim,  KbT  =  0.10 

stern  transom  trim,  KTrT  =  0.00 

initial  trim,  KT1T  =  0.00 

Trim  coefficient,  K_Tr  =  0.477 

par  Factors  and  constants 


Hul  1  i 

effect 

shallow,  Par 

_Hull_S 

= 

0. 

0133 

water 

depth  effect, 

par_ 

hT_S  = 

1. 

2431 

***  sinkage 

&  Trim 

Results  *** 

speed 

Fnh 

Par_ 

Par_ 

SMid 

Trim 

SbA 

SSA 

kts 

mps 

NO. 

Fnh 

hT 

(m/m) 

(m/m) 

(m) 

(m) 

5.00 

2.  57 

0.21 

0.05 

0.  91 

0.0010 

-0.0005 

0.  32 

0.19 

6.00 

3.09 

0.25 

0.07 

0.86 

0.0014 

-0.0007 

0.43 

0.26 

7.00 

3.  60 

0.29 

0.10 

0.  82 

0.0018 

-0.0008 

0.  56 

0.  35 

8.00 

4.12 

0.  34 

0.12 

0.77 

0.0023 

-0.0010 

0.71 

0.45 

9.00 

4.  63 

0.  38 

0.15 

0.  73 

0.0028 

-0.0012 

0.  86 

0.  57 

10.00 

5.14 

0.42 

0.18 

0.70 

0.0034 

-0. 0014 

1.04 

0.69 

11.00 

5.  66 

0.46 

0.  22 

0.  66 

0.0041 

-0.0015 

1.  22 

0.  84 

12.00 

6.17 

0.  50 

0.25 

0.63 

0.0048 

-0.0017 

1.43 

0.99 

13. 00 

6.  69 

0.  54 

0.  29 

0.  60 

0.0056 

-0.0019 

1.  65 

1.17 

14.00 

7.20 

0.  59 

0.  34 

0.  57 

0.0064 

-0.0021 

1.88 

1.  36 

15.00 

7.72 

0.63 

0.  39 

0.  55 

0.0074 

-0.0023 

2.14 

1.  56 

Figure  3b.  Example  generic  output  files  for  modified  Bunga  Saga  Emphat  Panamax  bulk  carrier 
and  unrestricted  channel  for  multiple  ship  speeds  AnkudMOut.out. 
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Figure  4.  Ankudinov  bow  and  stern  squat  predictions  for  modified  Bunga  Saga  Emphat  Panamax  bulk 

carrier  in  unrestricted  channel. 


Table  3.  ANKUDINOV4  and  ANKUDINOVM4  bow  and  stern  squat  predictions  for 

modified  Bunga  Saga  Emphat  Panamax  bulk  carrier. 

Channel  Type 

Location 

Unrestricted  U 

Restricted  R 

Canal  C 

5  knots 

10  knots 

5  knots 

10  knots 

5  knots 

10  knots 

Ankudinov  Predictions  (m) 

Bow 

0.32 

1.04 

0.30 

0.99 

0.35 

1.17 

Stern 

0.19 

0.69 

0.18 

0.68 

0.26 

0.94 

Ratio:  Ankudinov/PIANC 

Bow 

1.9 

1.5 

1.9 

1.5 

1.8 

1.5 

Stern 

2.4 

2.2 

2.0 

1.9 

2.0 

1.8 

The  modified  Bunga  Saga  Emphat  bulk  carrier  is  a  “blocky”  ship  (i.e.,  large  C«)  that  will  squat 
by  the  bow.  The  Ankudinov  predictions  are  on  the  conservative  side,  with  overpredictions 
approximately  2  times  (range  between  1.5  and  2.4)  those  of  other  PIANC  empirical  formulas.  In 
general,  the  Ankudinov  bow  predictions  decreased  as  the  ship  velocity  increased.  The  stern  pre- 
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dictions  also  decreased  as  speed  increased,  but  not  as  much  as  the  bow  predictions.  Bow  predic¬ 
tions  are  generally  better  (closer  to  PIANC  predictions)  than  the  stern  predictions,  which  is  good 
since  this  is  the  most  significant  squat  value  for  this  type  of  ship.  The  bow  predictions  ranged 
from  1.5  to  1.9  times  the  PIANC  values  compared  to  1.8  to  2.4  multipliers  for  the  stern  predic¬ 
tions.  Both  bow  and  stern  predictions  increased  for  the  C  channel  type  relative  to  the  U  channel, 
as  one  would  expect.  The  squat  for  the  R  channel  tended  to  decrease  slightly  or  stay  the  same 
relative  to  the  U  channel.  This  is  probably  due  to  the  fact  that  the  trench  height  was  only  25  per¬ 
cent  of  the  water  depth  and  the  side  slope  was  relatively  mild  for  such  a  wide  channel. 

Table  4  summarizes  the  minimum,  average,  and  maximum  ratios  for  the  three  channel  types  for 
the  range  of  speeds  from  5  to  15  knots.  For  bow  squat,  the  average  ratio  overpredictions  are  1.5 
for  all  three  channel  types.  Similarly  for  stern  squat,  the  average  squat  ratios  show  overpredic¬ 
tions  of  2.1,  1.8,  and  1.6  for  U,  R,  and  C  channels,  respectively.  Thus,  the  average  Ankudinov 
predictions  range  from  1.5  to  2  times  larger  than  the  PIANC  predictions  for  bow  and  stern  squat. 
In  summary,  this  comparison  with  PIANC  empirical  formulas  is  only  intended  to  show  the  rela¬ 
tive  values  of  the  Ankudinov  predictions  to  the  accepted  PIANC  standards.  Of  course,  the 
PIANC  values  are  also  only  predictions  and  not  measurements.  The  companion  Part  II  of  this 
CHETN  will  present  comparisons  and  validations  with  laboratory  and  field  measurements  for  a 
range  of  ship  and  channel  types. 


Table  4.  Statistical  ratios  of  Ankudinov  to  PIANC  bow  and  stern  squat  predictions  for 
modified  Bunga  Saga  Emphat  Panamax  bulk  carrier. 


Statistic 

Channel  Type 

Unrestricted  U 

Restricted  R 

Canal  C  j 

Bow 

Stern 

Bow 

Stern 

Bow 

Stern 

Minimum 

1.2 

1.8 

1.2 

1.4 

1.1 

0.7 

Average 

1.5 

2.1 

1.5 

1.8 

1.5 

1.6 

Maximum 

1.9 

2.4 

2.4 

2.0 

1.8 

2.0 

Notes: 

1 .  PIANC  Bow  squat  averages  based  on  12,  7,  and  6  values  for  U,  R,  and  C,  respectively. 

2.  PIANC  Stern  squat  averages  based  on  1  value  (Romisch)  for  all  three  channel  types. 

3.  Eleven  ship  speeds  from  5  to  15  knots  in  1  knot  increments  for  all  statistics. 


EXAMPLE  2:  A  second  example  is  included  of  the  Panamax  tanker  in  EM  1 1 10-2-1613  (Head¬ 
quarters,  U.S.  Army  Corps  of  Engineers  2006).  Table  5  lists  the  ship  and  channel  properties  for 
this  example  for  all  three  channel  types.  The  h/T  =1.1  for  this  example,  so  it  has  a  relatively 
shallow  UKC.  This  ship  will  squat  more  by  the  bow  since  Cb>  0.7.  Figure  5  compares 
the  Ankudinov  predictions  for  bow  squat  for  the  three  channel  types  with  those  of  PIANC’s 
Huuska  formula  from  the  EM  1110-2-1613  manual.  Table  6  lists  the  minimum,  average,  and 
maximum  ratios  of  Ankudinov  to  Huuska  bow  squat  for  the  three  channel  types.  These  compari¬ 
sons  with  the  PIANC  Huuska  predictions  are  closer  than  the  first  example  with  a  range  of  ratios 
of  0.9  to  2.0  for  all  three  channel  types.  The  Ankudinov  predictions  are  the  most  similar  to  the 
Huuska  predictions  for  the  canal  configuration  with  an  average  overprediction  of  1.1.  The  aver¬ 
age  overprediction  is  1.4  and  1.6  for  the  U  and  R  channels,  respectively. 
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Table  5.  Ship  and  channel  input  parameters  for  EM  111 0-2-1 61 3  Panamax  tanker. 

4p  (m) 

B(m) 

T  (m) 

Cs 

h  (m) 

Mm) 

W(m) 

n 

209.7 

32.3 

12.2 

0.85 

13.4 

7.32 

91.4 

3 

Notes: 

Tfp  =  Tap  =  T;  IProp  =  1;  IBow=  1;  IStern  =  0;  KpS  =  0.15;  KpT  =  0.15. 

Bow  Squat  for  EM1 1 10  Panamax  Tanker  -  U,  R  and  C 


AnkU 

AnkR 

AnkC 

Husll 


8 

V.,  knots 


10 


12 


Figure  5.  Ankudinov  bow  squat  predictions  for  EM  1 1 1 0-2-161 3  Panamax  tanker  in  U,  R,  and  C  channels. 


Table  6.  Statistical  ratios  of  Ankudinov  to  Huuska  bow  predictions  for  EM  110-2-1613 
Panamax  tanker. 

Statistic 

Channel  Type 

Unrestricted  U 

Restricted  R 

Canal  C 

Minimum 

1.1 

1.3 

0.9 

Average 

1.4 

1.6 

1.1 

Maximum 

1.8 

2.0 

1.3 
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SUMMARY:  This  technical  note  has  documented  the  Ankudinov  formulation  for  ship  squat  that 
was  used  in  the  ERDC  Ship/Tow  Simulator.  Two  FORTRAN  computer  programs  were  written: 
one  for  single  ship  speeds  and  one  for  multiple  ship  speeds.  The  theoretical  formulation  and  ship 
and  channel  input  and  output  parameters  are  described.  An  example  with  the  modified  Bunga 
Saga  Emphat  Panamax  bulk  carrier  that  was  used  in  CETN-I-63  and  CHETN-I-72  was  presented 
to  illustrate  the  input  and  output  parameters  and  comparisons  were  made  with  the  PIANC 
empirical  squat  formulas  for  bow  and  stem  squat  in  unrestricted,  restricted,  and  canal  channels. 
Finally,  another  example  for  the  EM  1 1 10-2-1613  Panamax  tanker  was  presented.  In  general,  the 
Ankudinov  formulas  tend  to  overpredict  bow  and  stem  squat  by  factors  of  two  and  larger  relative 
to  the  PIANC  empirical  formulas.  The  bow  predictions  were  better  (average  overprediction  of 
1.5)  than  the  stern  predictions,  which  is  good  since  both  example  ships  would  experience  maxi¬ 
mum  squat  by  the  bow.  The  PIANC  predictions  are  not  perfect,  but  are  generally  accepted  by  the 
port  and  harbor  community  as  reasonable  predictions  of  ship  squat.  In  Part  II  of  this  CHETN, 
further  comparisons  and  validations  with  laboratory  and  field  measurements  are  presented  and 
discussed. 

ADDITIONAL  INFORMATION:  This  CHETN  was  prepared  as  part  of  the  Deep  Draft  Naviga¬ 
tion  Research  Program,  Ship  Simulations  work  unit,  and  was  written  by  Dr.  Michael  J.  Briggs 
( Michael.J.Briggs@usace.army.mil ,  voice:  601-634-2005,  fax:  601-634-3433)  of  the  Coastal  and 
Hydraulics  Laboratory  (CHL),  U.S.  Army  Engineer  Research  and  Development  Center.  This 
technical  note  should  be  cited  as  follows: 

Briggs,  M.  J.  2009.  Ankudinov  ship  squat  predictions  -  Part  I:  Theory,  parame¬ 
ters,  and  FORTRAN  programs.  Coastal  and  Hydraulics  Engineering  Technical 
Note  ERDC/CHL  CHETN-IX-19,  Vicksburg,  MS:  U.S.  Army  Engineer  Research 
and  Development  Center,  http://chl.erdc.usace.army.mil/chetn/. 
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Appendix  A:  List  of  Symbols  in  FORTRAN  Programs  ANKUDINOV4  and 
ANKUDINOVM4 


Symbol 

Units 

Eq.  No. 

Description 

Equation 

Program 

BTr 

BTr 

... 

Stern  transom  width  used  in  KTrT 

F nh 

Fnh 

... 

6 

Depth  Froude  Number 

Kl 

KbT 

... 

15 

Bulbous  bow  factor 

KpS 

... 

3 

Propeller  sinkage  factor 

K 

KpT 

... 

14 

Propeller  trim  factor 

Krr 

KTr 

... 

12 

Trim  coefficient 

Kl 

KTrT 

... 

16 

Stern  transom  factor 

KTt, 

KT1T 

... 

17 

Initial  trim  effect  factor 

nTr 

nTr 

... 

13 

Trim  exponent 

P  Chi 

P_Ch1 

... 

8 

Channel  effect  parameter 

Pch2 

P_Ch2 

... 

18 

Channel  effect  trim  correction  parameter 

P Fnh 

Par_Fnh 

... 

5 

Ship  forward  speed  parameter 

PhT 

Par_hT_s 

... 

11 

Propeller  effect  in  shallow  water  on  trim  parameter 

Par_hT 

... 

7 

Water  depth  parameter 

Phu 

Par_Hull_s 

... 

4 

Ship  hull  parameter  for  shallow  water 

sb 

SbA 

m 

1 

Ship  squat  at  bow 

s„ 

Sh 

... 

9 

Channel  depth  factor  for  R  and  C  channels 

SMid 

m/m 

2 

Midpoint  ship  sinkage 

Ss 

SsA 

m 

1 

Ship  squat  at  stern 

Trim 

Trim 

m/m 

10 

Ship  trim 

Notes: 

1 .  Input  symbols  are  defined  in  Table  1 . 
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